Search for: All records

We quantify the cosmological spread of baryons relative to their initial neighbouring dark matter distribution using thousands of state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. We show that dark matter particles spread relative to their initial neighbouring distribution owing to chaotic gravitational dynamics on spatial scales comparable to their host dark matter halo. In contrast, gas in hydrodynamic simulations spreads much further from the initial neighbouring dark matter owing to feedback from supernovae (SNe) and active galactic nuclei (AGN). We show that large-scale baryon spread is very sensitive to model implementation details, with the fiducial simba model spreading ∼40 per cent of baryons >1 Mpc away compared to ∼10 per cent for the IllustrisTNG and astrid models. Increasing the efficiency of AGN-driven outflows greatly increases baryon spread while increasing the strength of SNe-driven winds can decrease spreading due to non-linear coupling of stellar and AGN feedback. We compare total matter power spectra between hydrodynamic and paired N-body simulations and demonstrate that the baryonic spread metric broadly captures the global impact of feedback on matter clustering over variations of cosmological and astrophysical parameters, initial conditions, and (to a lesser extent) galaxy formation models. Using symbolic regression, we find a function that reproduces the suppression of power by feedback as a function of wave number (k) and baryonic spread up to $k \sim 10\, h$ Mpc−1 in SIMBA while highlighting the challenge of developing models robust to variations in galaxy formation physics implementation.

 

An unprecedented array of new observational capabilities are starting to yield key constraints on models of the epoch of first light in the Universe. In this Letter we discuss the implications of the UV radiation background at cosmic dawn inferred by recent JWST observations for radio experiments aimed at detecting the redshifted 21 cm hyperfine transition of diffuse neutral hydrogen. Under the basic assumption that the 21 cm signal is activated by the Lyαphoton field produced by metal-poor stellar systems, we show that a detection at the low frequencies of the EDGES and SARAS3 experiments may be expected from a simple extrapolation of the declining UV luminosity density inferred atz≲ 14 from JWST early galaxy data. Accounting for an early radiation excess above the cosmic microwave background suggests a shallower or flat evolution to simultaneously reproduce low- and high-zcurrent UV luminosity density constraints, which cannot be entirely ruled out, given the large uncertainties from cosmic variance and the faint-end slope of the galaxy luminosity function at cosmic dawn. Our findings raise the intriguing possibility that a high star formation efficiency at early times may trigger the onset of intense Lyαemission at redshiftz≲ 20 and produce a cosmic 21 cm absorption signal 200 Myr after the Big Bang.

 

Abstract We train graph neural networks to perform field-level likelihood-free inference using galaxy catalogs from state-of-the-art hydrodynamic simulations of the CAMELS project. Our models are rotational, translational, and permutation invariant and do not impose any cut on scale. From galaxy catalogs that only contain 3D positions and radial velocities of ∼1000 galaxies in tiny ( 25 h − 1 Mpc ) 3 volumes our models can infer the value of Ω m with approximately 12% precision. More importantly, by testing the models on galaxy catalogs from thousands of hydrodynamic simulations, each having a different efficiency of supernova and active galactic nucleus feedback, run with five different codes and subgrid models—IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE—we find that our models are robust to changes in astrophysics, subgrid physics, and subhalo/galaxy finder. Furthermore, we test our models on 1024 simulations that cover a vast region in parameter space—variations in five cosmological and 23 astrophysical parameters—finding that the model extrapolates really well. Our results indicate that the key to building a robust model is the use of both galaxy positions and velocities, suggesting that the network has likely learned an underlying physical relation that does not depend on galaxy formation and is valid on scales larger than ∼10 h −1 kpc. 

We use two independent galaxy-formation simulations, flares, a cosmological hydrodynamical simulation, and shark, a semi-analytic model, to explore how well the JWST will be able to uncover the existence and parameters of the star-forming main sequence (SFS) at z = 5 → 10, i.e. shape, scatter, normalization. Using two independent simulations allows us to isolate predictions (e.g. stellar mass, star formation rate, SFR, luminosity functions) that are robust to or highly dependent on the implementation of the physics of galaxy formation. Both simulations predict that JWST can observe ≥70–90 per cent (for shark and flares, respectively) of galaxies up to z ∼ 10 (down to stellar masses of ${\approx}10^{8.3}\rm M_{\odot }$ and SFRs of ${\approx}10^{0.5}{\rm M}_{\odot }\,{\rm yr}^{-1}$) in modest integration times and given current proposed survey areas (e.g. the Web COSMOS 0.6 deg2) to accurately constrain the parameters of the SFS. Although both simulations predict qualitatively similar distributions of stellar mass and SFR. There are important quantitative differences, such as the abundance of massive, star-forming galaxies with flares predicting a higher abundance than shark; the early onset of quenching as a result of black hole growth in flares (at z ≈ 8), not seen in shark until much lower redshifts; and the implementation of synthetic photometry with flares predicting more JWST-detected galaxies (∼90 per cent) than shark (∼70 per cent) at z = 10. JWST observations will distinguish between these models, leading to a significant improvement upon our understanding of the formation of the very first galaxies.

 

ABSTRACT Matching the number counts of high-z submillimetre-selected galaxies (SMGs) has been a long-standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the submm emission from galaxies in the simba cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850 μm-selected sources. We find good agreement with the shape of the integrated 850 μm luminosity function, and the normalization is within 0.25 dex at >3 mJy, unprecedented for a fully cosmological hydrodynamic simulation, along with good agreement in the redshift distribution of bright SMGs. The agreement is driven primarily by simba’s good match to infrared measures of the star formation rate (SFR) function between z = 2 and 4 at high SFRs. Also important is the self-consistent on-the-fly dust model in simba, which predicts, on average, higher dust masses (by up to a factor of 2.5) compared to using a fixed dust-to-metals ratio of 0.3. We construct a light-cone to investigate the effect of far-field blending, and find that 52 per cent of sources are blends of multiple components, which makes a small contribution to the normalization of the bright end of the number counts. We provide new fits to the 850 μm luminosity as a function of SFR and dust mass. Our results demonstrate that solutions to the discrepancy between submm counts in simulations and observations, such as a top-heavy initial mass function, are unnecessary, and that submillimetre-bright phases are a natural consequence of massive galaxy evolution. 

The Baldwin, Philips, & Terlevich diagram of [Oiii]/Hβversus [Nii]/Hα(hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies atz∼ 2 reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from thesimbacosmological hydrodynamic galaxy formation simulation, using thecloudyphotoionization code to compute the nebular line luminosities from Hiiregions. We find that the observed shift toward higher [Oiii]/Hβand [Nii]/Hαvalues at high redshift arises from sample selection: when we consider only the most massive galaxiesM_*∼ 10^10–11M_⊙, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable toz∼ 0 observations. Even when accounting for samples-selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, Hiiregion densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies.

 

The SPT 0311–58 system atz= 6.900 is an extremely massive structure within the reionization epoch and offers a chance to understand the formation of galaxies at an extreme peak in the primordial density field. We present 70 mas Atacama Large Millimeter/submillimeter Array observations of the dust continuum and [Cii] 158μm emission in the central pair of galaxies and reach physical resolutions of ∼100–350 pc, among the most detailed views of any reionization-era system to date. The observations resolve the source into at least a dozen kiloparsec-size clumps. The global kinematics and high turbulent velocity dispersion within the galaxies present a striking contrast to recent claims of dynamically cold thin-disk kinematics in some dusty galaxies just 800 Myr later atz∼ 4. We speculate that both gravitational interactions and fragmentation from massive parent disks have likely played a role in the overall dynamics and formation of clumps in the system. Each clump individually is comparable in mass to other 6 <z< 8 galaxies identified in rest-UV/optical deep field surveys, but with star formation rates elevated by a factor of ~3-5. Internally, the clumps themselves bear close resemblance to greatly scaled-up versions of virialized cloud-scale structures identified in low-redshift galaxies. Our observations are qualitatively similar to the chaotic and clumpy assembly within massive halos seen in simulations of high-redshift galaxies.